suppressPackageStartupMessages(library(tidyverse))
devtools::load_all('~/Google Drive/My Drive/Scripts/R_packages/myUtilities/')
## ℹ Loading myUtilities
data_dir <- '/Volumes/Mitsu_NGS_3/METTL2A/'
wd <- "~/Google Drive/My Drive/Analysis/METTL2A/"
setwd(wd)
figdir <- paste0(wd, 'Figures/DRS_m3C_RNAs/Isoforms/Mapping_m3Csites_isoforms_ecdf_compare/')
tabledir <- paste0(wd, 'Tables/DRS_m3C_RNAs/Isoforms/')
theme_set(
theme_classic(base_size = 7) +
theme(legend.position = 'bottom')
)
read_signal_parameter_files <- function(
.glob,
.dir = '/Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/',
...) {
fs::dir_ls(.dir, glob = .glob) |>
read_tsv(
col_names = c('ref_pos', 'ref_kmer', 'dwell_time', 'intensity'),
col_types = 'dcdd',
id = 'filename', ...
) |>
separate(filename, sep = '_',
into = c(NA, NA, NA, NA, NA, 'sample_num', 'type', 'si', NA, 'transcript_id')) |>
mutate(transcript_id = str_remove(transcript_id, '.tsv'))
}
remove_outlier_from_eachsample <- function(df, outlier_thresh = .05) {
df |>
group_by(transcript_id, ref_pos, sample_num) |>
filter(intensity > quantile(intensity, probs = outlier_thresh / 2), na.rm = TRUE) |>
filter(intensity < quantile(intensity, probs = 1 - outlier_thresh / 2), na.rm = TRUE) |>
ungroup()
}
tsv_files <-
fs::dir_ls(
'/Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/',
glob = '*.tsv') |>
fs::path_filter(
regexp = 'ENST00000368716.9|ENST00000354332.8|ESPRESSO:chr1:1668:8'
)
tsv_files
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_7_siMETTL2A_I_N1_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_7_siMETTL2A_I_N1_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_8_siMETTL2A_I_N2_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_8_siMETTL2A_I_N2_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_9_siMETTL2A_I_N3_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221117_DrTaniue_9_siMETTL2A_I_N3_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_4_siMETTL2A_G_N1_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_4_siMETTL2A_G_N1_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_5_siMETTL2A_G_N2_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_5_siMETTL2A_G_N2_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_6_siMETTL2A_G_N3_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221122_DrTaniue_6_siMETTL2A_G_N3_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_1_Cont_D_N1_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_1_Cont_D_N1_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_2_Cont_D_N2_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_2_Cont_D_N2_ENST00000368716.9.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_3_Cont_D_N3_ENST00000354332.8.tsv
## /Volumes/Mitsu_NGS_2/METTL2A/Tables/Nanocompore/Signal_parameters/221128_DrTaniue_3_Cont_D_N3_ENST00000368716.9.tsv
parameters <-
tsv_files |>
map(read_signal_parameter_files) |>
reduce(bind_rows)
parameters
## # A tibble: 20,090,588 × 8
## sample_num type si transcript_id ref_pos ref_kmer dwell_time intensity
## <chr> <chr> <chr> <chr> <dbl> <chr> <dbl> <dbl>
## 1 7 siMETTL… I ENST00000354… 7 CTCTA 0.149 77.9
## 2 7 siMETTL… I ENST00000354… 8 TCTAC 0.0113 84.7
## 3 7 siMETTL… I ENST00000354… 9 CTACA 0.0150 83.5
## 4 7 siMETTL… I ENST00000354… 10 TACAA 0.00633 81.8
## 5 7 siMETTL… I ENST00000354… 11 ACAAC 0.00267 83.5
## 6 7 siMETTL… I ENST00000354… 12 CAACC 0.0120 89.2
## 7 7 siMETTL… I ENST00000354… 13 AACCC 0.0283 83.8
## 8 7 siMETTL… I ENST00000354… 14 ACCCT 0.00733 67.9
## 9 7 siMETTL… I ENST00000354… 15 CCCTC 0.00500 73.6
## 10 7 siMETTL… I ENST00000354… 16 CCTCT 0.00233 78.7
## # ℹ 20,090,578 more rows
S100A4_TCCTG <-
parameters |>
filter(ref_kmer == 'TCCTG') |>
filter(
(transcript_id == 'ENST00000368716.9' & ref_pos == 55) |
(transcript_id == 'ENST00000354332.8' & ref_pos == 100)
)
S100A4_TCCTG |>
export_tsv(outdir = tabledir, compression = 'gz')
##
## Exported to: ~/Google Drive/My Drive/Analysis/METTL2A/Tables/DRS_m3C_RNAs/Isoforms/S100A4_TCCTG_2024-07-29.tsv.gz
## # A tibble: 39,628 × 8
## sample_num type si transcript_id ref_pos ref_kmer dwell_time intensity
## <chr> <chr> <chr> <chr> <dbl> <chr> <dbl> <dbl>
## 1 7 siMETTL… I ENST00000354… 100 TCCTG 0.00467 78.7
## 2 7 siMETTL… I ENST00000354… 100 TCCTG 0.00267 73.7
## 3 7 siMETTL… I ENST00000354… 100 TCCTG 0.0200 72.8
## 4 7 siMETTL… I ENST00000354… 100 TCCTG 0.00300 73.1
## 5 7 siMETTL… I ENST00000354… 100 TCCTG 0.00967 74.8
## 6 7 siMETTL… I ENST00000354… 100 TCCTG 0.0110 74.1
## 7 7 siMETTL… I ENST00000354… 100 TCCTG 0.00633 76.3
## 8 7 siMETTL… I ENST00000354… 100 TCCTG 0.0110 77.1
## 9 7 siMETTL… I ENST00000354… 100 TCCTG 0.0163 71.9
## 10 7 siMETTL… I ENST00000354… 100 TCCTG 0.0100 74.8
## # ℹ 39,618 more rows
S100A4_TCCTG_intensity_ecdf <-
parameters |>
filter(ref_kmer == 'TCCTG') |>
filter(
(transcript_id == 'ENST00000368716.9' & ref_pos == 55) |
(transcript_id == 'ENST00000354332.8' & ref_pos == 100)
) |>
remove_outlier_from_eachsample() |>
ggplot(aes(x = intensity, colour = si, linetype = transcript_id, fill = ref_kmer)) +
scale_color_manual(values = c('#8C8C8C', '#37D9CC', '#A3A3F9')) +
stat_ecdf(lwd = 1.1)
S100A4_TCCTG_intensity_ecdf |>
ggsave_multiple_formats(
outdir = figdir, width = 6, height = 4, fontsize = 7
)
S100A4_ATCCT <-
parameters |>
filter(ref_kmer == 'ATCCT') |>
filter(
(transcript_id == 'ENST00000368716.9' & ref_pos == 54) |
(transcript_id == 'ENST00000354332.8' & ref_pos == 99)
)
S100A4_ATCCT |>
export_tsv(outdir = tabledir, compression = 'gz')
##
## Exported to: ~/Google Drive/My Drive/Analysis/METTL2A/Tables/DRS_m3C_RNAs/Isoforms/S100A4_ATCCT_2024-07-29.tsv.gz
## # A tibble: 39,634 × 8
## sample_num type si transcript_id ref_pos ref_kmer dwell_time intensity
## <chr> <chr> <chr> <chr> <dbl> <chr> <dbl> <dbl>
## 1 7 siMETTL… I ENST00000354… 99 ATCCT 0.0307 74.6
## 2 7 siMETTL… I ENST00000354… 99 ATCCT 0.00767 71.1
## 3 7 siMETTL… I ENST00000354… 99 ATCCT 0.00733 72.3
## 4 7 siMETTL… I ENST00000354… 99 ATCCT 0.0123 72.5
## 5 7 siMETTL… I ENST00000354… 99 ATCCT 0.00700 70.9
## 6 7 siMETTL… I ENST00000354… 99 ATCCT 0.0107 73.9
## 7 7 siMETTL… I ENST00000354… 99 ATCCT 0.00300 73.3
## 8 7 siMETTL… I ENST00000354… 99 ATCCT 0.00433 76.5
## 9 7 siMETTL… I ENST00000354… 99 ATCCT 0.00267 71.9
## 10 7 siMETTL… I ENST00000354… 99 ATCCT 0.00267 72.0
## # ℹ 39,624 more rows
S100A4_ATCCT_intensity_ecdf <-
parameters |>
filter(ref_kmer == 'ATCCT') |>
filter(
(transcript_id == 'ENST00000368716.9' & ref_pos == 54) |
(transcript_id == 'ENST00000354332.8' & ref_pos == 99)
) |>
remove_outlier_from_eachsample() |>
ggplot(aes(x = intensity, colour = si, linetype = transcript_id, fill = ref_kmer)) +
scale_color_manual(values = c('#8C8C8C', '#37D9CC', '#A3A3F9')) +
stat_ecdf(lwd = 1.1)
S100A4_ATCCT_intensity_ecdf |>
ggsave_multiple_formats(
outdir = figdir, width = 6, height = 4, fontsize = 7
)
Density
S100A4_TCCTG_intensity_density <-
parameters |>
filter(ref_kmer == 'TCCTG') |>
filter(
(transcript_id == 'ENST00000368716.9' & ref_pos == 55) |
(transcript_id == 'ENST00000354332.8' & ref_pos == 100)
) |>
remove_outlier_from_eachsample(outlier_thresh = .1) |>
ggplot(aes(x = intensity, fill = si, linetype = transcript_id, colour = si)) +
scale_fill_manual(values = c('#8C8C8C', '#37D9CC', '#A3A3F9')) +
scale_color_manual(values = c('#8C8C8C', '#37D9CC', '#A3A3F9')) +
geom_density(alpha = 1/2) +
facet_wrap( ~ transcript_id, ncol = 1)
S100A4_TCCTG_intensity_density |>
ggsave_multiple_formats(
outdir = figdir, width = 6, height = 4, fontsize = 7
)
## Warning in grid.Call.graphics(C_polygon, x$x, x$y, index): semi-transparency is
## not supported on this device: reported only once per page
parameters |>
filter(ref_kmer == 'TCCTG') |>
filter(
#(transcript_id == 'ENST00000368716.9' & ref_pos == 55)
(transcript_id == 'ENST00000354332.8' & ref_pos == 100)
) |>
remove_outlier_from_eachsample() |>
ggplot(aes(
x = intensity, colour = si, linetype = sample_num, fill = ref_kmer)) +
scale_color_manual(values = c('#8C8C8C', '#37D9CC', '#A3A3F9')) +
stat_ecdf(lwd = 1.1)
sessioninfo::session_info()
## ─ Session info ───────────────────────────────────────────────────────────────
## setting value
## version R version 4.3.1 (2023-06-16)
## os macOS Sonoma 14.5
## system aarch64, darwin20
## ui X11
## language (EN)
## collate en_US.UTF-8
## ctype en_US.UTF-8
## tz Asia/Tokyo
## date 2024-07-29
## pandoc 3.1.1 @ /Applications/RStudio.app/Contents/Resources/app/quarto/bin/tools/ (via rmarkdown)
##
## ─ Packages ───────────────────────────────────────────────────────────────────
## ! package * version date (UTC) lib source
## abind 1.4-5 2016-07-21 [1] CRAN (R 4.3.0)
## backports 1.4.1 2021-12-13 [1] CRAN (R 4.3.0)
## bit 4.0.5 2022-11-15 [1] CRAN (R 4.3.0)
## bit64 4.0.5 2020-08-30 [1] CRAN (R 4.3.0)
## broom 1.0.5 2023-06-09 [1] CRAN (R 4.3.0)
## bslib 0.5.1 2023-08-11 [1] CRAN (R 4.3.0)
## cachem 1.0.8 2023-05-01 [1] CRAN (R 4.3.0)
## callr 3.7.3 2022-11-02 [1] CRAN (R 4.3.0)
## car 3.1-2 2023-03-30 [1] CRAN (R 4.3.0)
## carData 3.0-5 2022-01-06 [1] CRAN (R 4.3.0)
## class 7.3-22 2023-05-03 [1] CRAN (R 4.3.1)
## cli 3.6.1 2023-03-23 [1] CRAN (R 4.3.0)
## codetools 0.2-19 2023-02-01 [1] CRAN (R 4.3.1)
## colorspace 2.1-0 2023-01-23 [1] CRAN (R 4.3.0)
## crayon 1.5.2 2022-09-29 [1] CRAN (R 4.3.0)
## data.table 1.14.8 2023-02-17 [1] CRAN (R 4.3.0)
## desc 1.4.2 2022-09-08 [1] CRAN (R 4.3.0)
## devtools 2.4.5 2022-10-11 [1] CRAN (R 4.3.0)
## dials 1.2.0 2023-04-03 [1] CRAN (R 4.3.0)
## DiceDesign 1.9 2021-02-13 [1] CRAN (R 4.3.0)
## digest 0.6.33 2023-07-07 [1] CRAN (R 4.3.0)
## dplyr * 1.1.2 2023-04-20 [1] CRAN (R 4.3.0)
## ellipsis 0.3.2 2021-04-29 [1] CRAN (R 4.3.0)
## evaluate 0.21 2023-05-05 [1] CRAN (R 4.3.0)
## fansi 1.0.4 2023-01-22 [1] CRAN (R 4.3.0)
## farver 2.1.1 2022-07-06 [1] CRAN (R 4.3.0)
## fastmap 1.1.1 2023-02-24 [1] CRAN (R 4.3.0)
## forcats * 1.0.0 2023-01-29 [1] CRAN (R 4.3.0)
## foreach 1.5.2 2022-02-02 [1] CRAN (R 4.3.0)
## fs 1.6.3 2023-07-20 [1] CRAN (R 4.3.0)
## furrr 0.3.1 2022-08-15 [1] CRAN (R 4.3.0)
## future 1.33.0 2023-07-01 [1] CRAN (R 4.3.0)
## future.apply 1.11.0 2023-05-21 [1] CRAN (R 4.3.0)
## generics 0.1.3 2022-07-05 [1] CRAN (R 4.3.0)
## ggforce 0.4.1 2022-10-04 [1] CRAN (R 4.3.0)
## ggplot2 * 3.4.3 2023-08-14 [1] CRAN (R 4.3.0)
## ggpubr 0.6.0 2023-02-10 [1] CRAN (R 4.3.0)
## ggrepel 0.9.3 2023-02-03 [1] CRAN (R 4.3.0)
## ggsignif 0.6.4 2022-10-13 [1] CRAN (R 4.3.0)
## globals 0.16.2 2022-11-21 [1] CRAN (R 4.3.0)
## glue 1.6.2 2022-02-24 [1] CRAN (R 4.3.0)
## gower 1.0.1 2022-12-22 [1] CRAN (R 4.3.0)
## GPfit 1.0-8 2019-02-08 [1] CRAN (R 4.3.0)
## gprofiler2 0.2.3 2024-02-23 [1] CRAN (R 4.3.1)
## gtable 0.3.3 2023-03-21 [1] CRAN (R 4.3.0)
## hardhat 1.3.0 2023-03-30 [1] CRAN (R 4.3.0)
## highr 0.10 2022-12-22 [1] CRAN (R 4.3.0)
## hms 1.1.3 2023-03-21 [1] CRAN (R 4.3.0)
## htmltools 0.5.6 2023-08-10 [1] CRAN (R 4.3.0)
## htmlwidgets 1.6.2 2023-03-17 [1] CRAN (R 4.3.0)
## httpuv 1.6.11 2023-05-11 [1] CRAN (R 4.3.0)
## httr 1.4.7 2023-08-15 [1] CRAN (R 4.3.0)
## ipred 0.9-14 2023-03-09 [1] CRAN (R 4.3.0)
## iterators 1.0.14 2022-02-05 [1] CRAN (R 4.3.0)
## jquerylib 0.1.4 2021-04-26 [1] CRAN (R 4.3.0)
## jsonlite 1.8.7 2023-06-29 [1] CRAN (R 4.3.0)
## knitr 1.43 2023-05-25 [1] CRAN (R 4.3.0)
## labeling 0.4.2 2020-10-20 [1] CRAN (R 4.3.0)
## later 1.3.1 2023-05-02 [1] CRAN (R 4.3.0)
## lattice 0.21-8 2023-04-05 [1] CRAN (R 4.3.1)
## lava 1.7.2.1 2023-02-27 [1] CRAN (R 4.3.0)
## lazyeval 0.2.2 2019-03-15 [1] CRAN (R 4.3.0)
## lhs 1.1.6 2022-12-17 [1] CRAN (R 4.3.0)
## lifecycle 1.0.3 2022-10-07 [1] CRAN (R 4.3.0)
## listenv 0.9.0 2022-12-16 [1] CRAN (R 4.3.0)
## lubridate * 1.9.2 2023-02-10 [1] CRAN (R 4.3.0)
## magrittr 2.0.3 2022-03-30 [1] CRAN (R 4.3.0)
## MASS 7.3-60 2023-05-04 [1] CRAN (R 4.3.1)
## Matrix 1.6-1 2023-08-14 [1] CRAN (R 4.3.0)
## memoise 2.0.1 2021-11-26 [1] CRAN (R 4.3.0)
## mime 0.12 2021-09-28 [1] CRAN (R 4.3.0)
## miniUI 0.1.1.1 2018-05-18 [1] CRAN (R 4.3.0)
## munsell 0.5.0 2018-06-12 [1] CRAN (R 4.3.0)
## R myUtilities * 0.0.0.9000 <NA> [?] <NA>
## nnet 7.3-19 2023-05-03 [1] CRAN (R 4.3.1)
## parallelly 1.36.0 2023-05-26 [1] CRAN (R 4.3.0)
## parsnip 1.1.1 2023-08-17 [1] CRAN (R 4.3.0)
## pillar 1.9.0 2023-03-22 [1] CRAN (R 4.3.0)
## pkgbuild 1.4.2 2023-06-26 [1] CRAN (R 4.3.0)
## pkgconfig 2.0.3 2019-09-22 [1] CRAN (R 4.3.0)
## pkgload 1.3.2.1 2023-07-08 [1] CRAN (R 4.3.0)
## plotly 4.10.4 2024-01-13 [1] CRAN (R 4.3.1)
## polyclip 1.10-4 2022-10-20 [1] CRAN (R 4.3.0)
## prettyunits 1.1.1 2020-01-24 [1] CRAN (R 4.3.0)
## processx 3.8.2 2023-06-30 [1] CRAN (R 4.3.0)
## prodlim 2023.08.28 2023-08-28 [1] CRAN (R 4.3.0)
## profvis 0.3.8 2023-05-02 [1] CRAN (R 4.3.0)
## promises 1.2.1 2023-08-10 [1] CRAN (R 4.3.0)
## ps 1.7.5 2023-04-18 [1] CRAN (R 4.3.0)
## purrr * 1.0.2 2023-08-10 [1] CRAN (R 4.3.0)
## R6 2.5.1 2021-08-19 [1] CRAN (R 4.3.0)
## ragg 1.2.5 2023-01-12 [1] CRAN (R 4.3.0)
## Rcpp 1.0.11 2023-07-06 [1] CRAN (R 4.3.0)
## readr * 2.1.4 2023-02-10 [1] CRAN (R 4.3.0)
## recipes 1.0.8 2023-08-25 [1] CRAN (R 4.3.0)
## remotes 2.4.2.1 2023-07-18 [1] CRAN (R 4.3.0)
## rlang 1.1.1 2023-04-28 [1] CRAN (R 4.3.0)
## rmarkdown 2.24 2023-08-14 [1] CRAN (R 4.3.0)
## rpart 4.1.19 2022-10-21 [1] CRAN (R 4.3.1)
## rprojroot 2.0.3 2022-04-02 [1] CRAN (R 4.3.0)
## rsample 1.2.0 2023-08-23 [1] CRAN (R 4.3.0)
## rstatix 0.7.2 2023-02-01 [1] CRAN (R 4.3.0)
## rstudioapi 0.15.0 2023-07-07 [1] CRAN (R 4.3.0)
## sass 0.4.7 2023-07-15 [1] CRAN (R 4.3.0)
## scales 1.2.1 2022-08-20 [1] CRAN (R 4.3.0)
## sessioninfo 1.2.2 2021-12-06 [1] CRAN (R 4.3.0)
## shiny 1.7.5 2023-08-12 [1] CRAN (R 4.3.0)
## stringi 1.7.12 2023-01-11 [1] CRAN (R 4.3.0)
## stringr * 1.5.0 2022-12-02 [1] CRAN (R 4.3.0)
## survival 3.5-5 2023-03-12 [1] CRAN (R 4.3.1)
## svglite 2.1.1 2023-01-10 [1] CRAN (R 4.3.0)
## systemfonts 1.0.4 2022-02-11 [1] CRAN (R 4.3.0)
## textshaping 0.3.6 2021-10-13 [1] CRAN (R 4.3.0)
## tibble * 3.2.1 2023-03-20 [1] CRAN (R 4.3.0)
## tidyr * 1.3.0 2023-01-24 [1] CRAN (R 4.3.0)
## tidyselect 1.2.0 2022-10-10 [1] CRAN (R 4.3.0)
## tidyverse * 2.0.0 2023-02-22 [1] CRAN (R 4.3.0)
## timechange 0.2.0 2023-01-11 [1] CRAN (R 4.3.0)
## timeDate 4022.108 2023-01-07 [1] CRAN (R 4.3.0)
## tune 1.1.2 2023-08-23 [1] CRAN (R 4.3.0)
## tweenr 2.0.2 2022-09-06 [1] CRAN (R 4.3.0)
## tzdb 0.4.0 2023-05-12 [1] CRAN (R 4.3.0)
## urlchecker 1.0.1 2021-11-30 [1] CRAN (R 4.3.0)
## usethis 2.2.2 2023-07-06 [1] CRAN (R 4.3.0)
## utf8 1.2.3 2023-01-31 [1] CRAN (R 4.3.0)
## vctrs 0.6.3 2023-06-14 [1] CRAN (R 4.3.0)
## viridisLite 0.4.2 2023-05-02 [1] CRAN (R 4.3.0)
## vroom 1.6.3 2023-04-28 [1] CRAN (R 4.3.0)
## withr 2.5.0 2022-03-03 [1] CRAN (R 4.3.0)
## workflows 1.1.3 2023-02-22 [1] CRAN (R 4.3.0)
## xfun 0.40 2023-08-09 [1] CRAN (R 4.3.0)
## xtable 1.8-4 2019-04-21 [1] CRAN (R 4.3.0)
## yaml 2.3.7 2023-01-23 [1] CRAN (R 4.3.0)
## yardstick 1.2.0 2023-04-21 [1] CRAN (R 4.3.0)
##
## [1] /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/library
##
## R ── Package was removed from disk.
##
## ──────────────────────────────────────────────────────────────────────────────